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MEASUREMENT OF EFFECTIVE DELAYED NEUTRON FRACTION 
FOR NASA ZERO POWER REACTOR I 
by Donald Alger, Wendell Mayo, and Robert Mueller 
Lewis Research Center 

SUMMARY 

The effective delayed neutron fraction /3 e ^ has been experimentally determined for 
each of three uranyl fluoride - water solution reactors in the NASA Zero Power Reactor I 
facility. The hydrogen to uranium 235 atom ratios R for the solutions in these measure- 
ments are 190, 473, and 565. The value of /3 e ^ has been measured by the boron substi- 
tution method. Values for ( 3 ^ obtained are 0. 0090±0. 0006, 0. 0086±0. 0005, and 
0. 0082±0. 0004 for the respective atom ratios. Calculated values of /3 e ^ are in reason- 
able agreement with the experimental results. 


INTRODUCTION 

The inhour method (ref. 1) is a kinetic technique that is commonly used to measure 
the reactivity associated with perturbations to a delayed critical reactor. The reactivity 
is measured in terms of the effective delayed neutron fraction This effective de- 

layed neutron fraction is usually chosen as the unit of reactivity and called the dollar 
(ref. 2). 

The effective delayed neutron fraction for a given reactor differs from the de- 

layed neutron fraction /3 (ref. 14) because the delayed neutrons are emitted with a lower 
average energy (~0. 5 MeV) than the prompt neutrons (~2 MeV). Delayed neutrons there- 
fore have a smaller leakage probability, and their effectiveness in causing new fissions in 
a thermal reactor core is greater than the prompt neutrons. Effects on due to ura- 

nium 238 fast fissions do not occur since nearly all of the delayed neutrons are emitted 
below the uranium 238 fission energy threshold (~0. 8 MeV). Therefore, interaction of 
delayed neutrons with fast fissions in uranium 238 is a very small effect in uranyl fluoride 
solution reactors that are essentially thermal reactors. Values of or the ratios of 



/3 e ff to jS(sy), are functions of the composition and geometry of the reactors involved 
(refs. 3 and 4). 

The NASA Zero Power Reactor I (ZPR-I) is a uranyl fluoride - water (UOgFg-HgO) 
system with criticality achieved by varying solution height. A wide range of uranium 235 
loadings is possible depending on the solution density. A measure of the loading is the 
hydrogen to uranium 235 atom ratio R. Criticality in the unreflected system may be 
achieved with less than 3 kilograms of uranium 235 with R values ranging from approxi- 
mately 180 to 570. 

There were two objectives of this study, one of which was to determine /3 e ^ for 
several uranyl fluoride - water solution concentrations. (A calibration curve relating the 
primary measurement, or the change in fuel solution height, to absolute reactivity may 
then be constructed. ) The other objective was to test calculational techniques that may be 
used to extend these results to other reactor configurations. 

Gwin and Weinberg (ref. 5) determined <3^ for several uranyl fluoride - water solu- 
tion reactors using experimentally determined bucklings. The measurement of was 

not the primary objective of that study and the results were approximate. The method re- 
quired reactors of various solid configurations (spheres, cylinders, etc.). Ruane, et al. 
(ref. 6), Kaplan and Henry (ref. 7), and Perez-Belles, Kington, and DeSaussure (ref. 8) 
used the boron substitution technique for the measurement of (3 q ^ in their respective re- 
actors. The measurements reported herein for the ZPR-I reactor follow the methods of 
reference 8. The calculational method employed follows that reported in reference 7. hi 
both cases, some modification and extension of the theory is necessary for application to 
the ZPR-I reactor. 

The following section presents the theory of the experiment. Next follow the experi- 
mental procedure, the description of the equipment, and the measurements and results, 
followed by the theoretical calculations for /3 eff . 

THEORY OF THE EXPERIMENT 

The theory of the experiment is based on two expressions for the reactivity due to 
perturbations to a critical reactor, the inhour equation, and the time -independent pertur- 
bation equation. Application of these two equations yields an expression for j3^ which, 
in many reactors, is well suited to experimental evaluation. The results have been ap- 
plied, for example, to the Bulk Shielding Reactor-I at the Oak Ridge National Laboratory 
(ref. 8). The notation used in reference 8 is used here when applicable. The ZPR-I 
measurements for /3 g ^ depend on modification and extension of the results given in ref- 
erence 8. The reasons for this and the problems presented are discussed as the theory 
is developed. 
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It is assumed that a nonfissionable material exists whose scattering and absorption 
cross sections at every energy are identical with those of the fuel. If a small amount of 
this material is substituted for an equal volume of fuel in a critical reactor, the per- 
turbed reactor will have a reactivity given by the inhour equation 


6p 


T i 




(1 + *jT.) 


( 1 ) 


Symbols are defined in appendix A. 

It is assumed that the substitution in a small region of volume results in a long 
stable reactor period Tp which is greater than 200 seconds for these experiments. The 
mean prompt neutron lifetime i for the ZPR-I varies considerably with fuel solution 
concentration, but it is of the order of 40 microseconds for the range of concentrations 
studied. Thus, the term £/Tj in equation (1) may be neglected for these measurements. 

The same reactivity for a perturbed reactor is given by the time -independent pertur- 
bation equation 


6p = 


/ dr/dE6 !vS f (r, E)®(r, E)P*(r) 


d I E)®(r, E)P*(r) - J dr JdEvZ f (r, E)$(r, E)P*(r) 


( 2 ) 


where the fission spectrum weighted neutron importance is 


P*(r) 



E’) 


In principle, only a single substitution is necessary since the measured quantity is 
the reactor period; hence, if equation (2) can be evaluated with sufficient accuracy, 
can be found from equations (1) and (2). Equation (2) may be evaluated in several ways, 
depending on the assumptions made. Henry (ref. 9) has shown that two static multigroup 
spatial solutions in which the geometry is accurately represented can be used to obtain 
6p T fully in agreement with equation (2). In performing a substitution experiment for 

the effective delayed neutron fraction in a solid fueled reactor, Kaplan and Henry (ref. 7) 
chose to sum the results obtained experimentally for successive substitutions in different 
fueled regions and to sum the results obtained by evaluation of equation (2) using four 
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energy group spatial flux and adjoint calculated solutions to get the best value for 
that is, 



If equation (1) is used, it follows that 



where i represents the i** 1 substitution and 


sign/6p T \ 

' i'calc 


is negative when 6p T is 


negative and positive when 6 p is positive. 

If, however, successive substitutions are made over the entire fueled region of the 
core, or at least over symmetric regions, equation (3) reduces to an expression involving 
only experimentally determined periods. To show this, equation (2) is summed over all 
regions Tj that make up the total core volume V. A period measurement is made for 
each substitution and care is taken not to overlap regions. The substitutions must result 
in a small but well defined reactivity change. Equation (2) is now written as 


Z 6 v 


J dr JdE6 vY f (r, E)<f>(r, E)P*(r) 


dr JdEvY f (r, E)$(r, E)P*(r) - dr J'^EvY^r, E)$(r, E)P*(r)^ 


( 4 ) 


where the second term in the denominator of equation (4) is an average for all substitu- 
tions of that quantity. Since, for these experiments, the second term in the denominator 
of equation (4) is small compared to the first term, it may be neglected. For these ex- 
periments, it follows that 


Si'Sj, = -vY^ 
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mi in hi 


■■■ I III 


and therefore 



(5) 


Hence, equation (3) may be written as 



( 6 ) 


A correction to equations (5) and (6) usually must be made to account for the mis- 
match in cross sections between the fuel and the substitute. If the materials for the sub- 
stitute are chosen properly, then the only mismatch of any consequence is due to the ab- 
sorption cross section. The correction may be obtained by using the absorption term 
from the time-independent perturbation equation as 


6A = 


J dr ydE62 a (r, E)$(r, E)$*(r, E) 
/ dr JdEvZ t (r, E)3?(r, E)P*(r) 


(7) 


where the second term in the denominator is omitted. The correction is usually small 
and, for these experiments, did not exceed 0. 2 percent. Equation (6) is now rewritten to 
include this corrected 6A as 



( 8 ) 


Because the ZPR-I utilizes the uranyl fluoride fuel solution, it is impractical to at- 
tempt substitutions that map out the entire fuel volume, or even a symmetric segment of 
the fuel. Use is made of the fact noted in equation (2) that the reactivity is proportional 
to the product of the power and fission spectrum weighted adjoint; that is. 
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(9) 


6p T cc J' dr JdE i>S f (r, E)$(r, E)P*(r) = f P(r)P*(r)dr 


where 


P(r) = fc iEvS f (r, E)$(r, E) 

Several substitutions of small, but different, volumes at the same position in the core 

allow the experimental specific reactivity p to be determined, and equation (8) becomes 

s 


0 , 


1 - 6A 


eff 


— / P(r)P*(r)dr 

R Jv ~ ~ ~ 


(10) 


where the 


P(r) and P*(r) are normalized such that 


Y' .. Vf . f P(r)P*(r)dr 

Z^i +x jTi P \ “ “ " 


( 11 ) 


If the extrapolation distances at the boundaries of an unreflected reactor are assumed 
to be energy independent, then P(r) and P*(r) have the same spatial distribution; <f»(r) 
and $*(r) also have the same spatial distributions, where 


and 


-/ 


$(r) = / dE$(r, E) 


<f>*(r) = J dE$*(r, E) 


Equation (10) is rewritten as 


^eff - - 1-SA < 12 > 

— /$( r)$*(r)dr 

0 J ~ ~ ~ 


where the integral in the denominator is normalized as in equation (11). This expression 
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is used for the determination of . The experimental and computational procedures 
are described in the following sections. 


EXPERIMENTAL PROCEDURE 
Reactor Description 


The ZPR-I is an unreflected (bare) cylindrical reactor utilizing a uranyl fluoride - 
water solution as fuel and moderator. It is, therefore, a homogeneous thermal reactor. 
Criticality, for a given fuel concentration, is achieved by varying the fuel solution height. 
The height is measured by a micrometer lead screw with a platinum electrical continuity 
probe attached. The uranium is enriched to 93. 2 percent uranium 235. The aluminum 
tank has a 30. 48-centimeter inside diameter with a 0. 3 175 -centimeter -thick wall. The 
critical height varies between approximately 23 and 76 centimeters depending on the fuel 
solution concentration. 

Direct measurement of the stable reactor period that would result if a neutron ab- 
sorber were inserted or removed from a ZPR-I core is not possible, since no means are 
available for the remote movement of an experiment into or out of an already critical 
ZPR-I core. Therefore, the reactivity associated with a neutron absorber placed into a 
ZPR-I core must be obtained indirectly by measuring the critical height of an unperturbed 
core, shutting the reactor down by draining the solution, and then positioning the absorber 
in the core and measuring the new critical height of the solution. The difference in criti- 
cal height measurements must then be related to reactivity by a previously obtained cali- 
bration curve. A ZPR-I calibration curve, and a discussion of its determination, is pre- 
sented in appendix B. 

The ZPR-I base cores chosen for these experiments had critical heights of 24. 892 
centimeters (9. 800 in. ), 45. 720 centimeters (18. 000 in. ), and 67. 310 centimeters 


TABLE I. - FUEL SOLUTION CONSTANTS OF ZERO POWER REACTOR I BASE CORES 


Base core 

Uranyl 

Uranyl fluoride 

- water atom concentration, 

Hydrogen to 

Macroscopic 

critical 

fluoride - 



N, 



uranium 235 

neutron 

height 
at 20° C, 
cm 

water 

solution 

density 



atoms/cm 6 



atom ratio, 

absorption 
cross section, 
2.(0.0253 eV), 

Hydrogen 

Uranium 235 

Uranium 238 

Oxygen 

Fluorine 

R 


at 20° C, 







cm" 1 


gm/cm 3 







(a) 

24.892 

1.1672 

6. 5505X10 22 

3.4463X10 20 

0.2511X10 20 

3. 3492X10 22 

7. 3947X10 20 

190.1 

0.2609 

45.720 

1.0671 

6.6263 

1.4005 

. 1021 

3.3432 

3.0052 

473.1 

.1191 

67.310 

1.0559 

6.6346 

1. 1726 

.08545 

3.3424 

2. 5161 

565.8 

.1033 


a o A (0. 0253 eV) from table U. 
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TABLE n. - MICROSCOPIC neutron ( 26 500 ) Fuel solution constants for these base 

absorption cross sections cores are given in table I (p. 7). The macroscopic 

neutron absorption cross sections S A (0. 0253 eV) are 
calculated from the atom concentrations of table I, 
and the microscopic neutron absorption cross sections 
ct a ( 0. 0253 eV) in table II. 

The critical height of a solution reactor decreases 
slightly from day to day because of the evaporation of 
water from the fuel solution. It is not feasible, how- 
ever, to add water each day to make up for the evapo- 
ration loss. As a result, the unperturbed critical 
solution height during a series of measurements does 
not remain constant. Evaporation of water from the 
uranyl fluoride - water fuel solution is, however, essentially a uniform perturbation of 
the solution S A (0. 0253 eV) with a certain reactivity worth relative to some chosen base 
critical height. Each of the experiments can be corrected to the base core height by use 
of the calibration curve in appendix B. 


Element 

Microscopic neutron 
absorption 
cross section, 
a A (0.0253 eV), 
b 

Hydrogen 

0.33 

Uranium 235 

694.0 

Uranium 238 

2.71 

Oxygen 

20X10“ 5 

Fluorine 

.01 

Boron (natural) 

762.0 


Boron -for -Fuel Substitution 


Pure analytical reagent grade boric acid crystals dissolved in water were used as a 
substitute for fuel. The boron content of the boric acid solutions was specified by equat- 
ing the macroscopic thermal neutron absorption cross section ( S A (0. 0253 eV)) of a par- 
ticular uranyl fluoride - water solution, in which the substitution was to be made, to that 
of the boric acid solution. The boric acid solutions were prepared at Lewis. 

Five thin-walled plastic spheres were fabricated. Volumes of the spheres are shown 
in table in. A 0. 1588-centimeter-diameter (1/16-in. -diam) filling hole was tapped in the 

surface of each ball and sealed with a threaded 
plastic screw and a neoprene O-ring. Vol- 
umes were determined by filling these spheres 
with 20° C demineralized water, weighing, 
and then calculating each volume from the 
mass of water and known absolute density of 
water at 20° C (0. 998203 g/cm^). In deter - 
minining the volumes, care was taken to re- 
move all air bubbles. 

For the experiments, the spheres were 
filled with either boric acid solution or fuel 


TABLE m. - VOLUMES OF SPHERICAL PLASTIC 
CONTAINERS USED FOR BORON-FOR-FUEL 
SUBSTITUTION MEASUREMENTS 


Plastic ball number 

Ball volume, 


cm ^ 

1 

28. 3000 

2 

19. 3249 

3 

9.4124 

4 

6.3727 

5 

4.6477 
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II 





solution and positioned in the core by a stainless -steel wire holder. A sketch of the 
ZPR-I with a plastic sphere held in the core center is shown in figure 1 (p. 9). 

For each base core, three concentrations of boric acid solution were prepared that 
would bracket the solution 2^(0. 0253 eV). This is necessary since the 2^(0.0253 eV) 
of the uranyl fluoride - water solution varies slightly with time because of evaporation of 
water from the solution. Prior to a substitution experiment, a sample of ZPR-I solution 
was removed and its density measured at 20° C. Uranium 235 and other constituent atom 
densities N, in atoms per cubic centimeter of solution, were determined from previously 
obtained experimental data relating fuel solution density at 20° C and uranium atom con- 
tent as determined by a chemical analysis (ref . 10). Atom densities, microscopic and 
macroscopic absorption cross sections for the fuel, and boric acid solutions are given in 
tables I, n, and IV. 

Measurement of the boron-for-fuel substitution reactivity worth for each base core 
involved first filling the five plastic balls (see table m, p. 8) with uranyl fluoride - water 
solution of the base core. No attempt was made to remove completely all air bubbles 
from the balls when they were filled with fuel. Consequently, a small air bubble re- 
mained. This filled volume was calculated from the mass of fuel in a ball and the mea- 
sured fuel density. Air bubble volumes equalled the difference between the fuel volume 
in the ball and the volumes in table HI. After the measurement of the reactivity worths 

TABLE IV. - DATA FOR BORIC ACID SOLUTION CONCENTRATIONS 
USED AS BORON-FOR-FUEL SUBSTITUTIONS IN 
ZERO POWER REACTOR I BASE CORES 


Base core solution 

Concentration 
of boric acid 

Boric acid solution 





Critical 

Macroscopic 

solution, 

Density, 

Macroscopic 

height 
at 20° C, 
cm 

neutron absorption 
cross section 
at 0.0253 eV, 

S A (0.0253 eV), 

A -1 
cm 

q 

g solute/cm 

gm/cm** 

neutron absorption 
cross section 
at 0.0253 eV, 

S A (0.0253 eV), 

A -1 
cm 

24. 892 

0.2609 

0.028987 

1.00888 

0.2371 


.2609 

.031938 

1.00993 

.2590 


.2609 

.035008 

1.01111 

.2817 

45.720 

0. 1191 

0.011579 

1.00246 

0. 1079 


. 1191 

.013602 

1.00317 

.1229 


.1191 

.016040 

1.00410 

.1410 

67.310 

0. 1033 

0.01001 

1.00189 

0. 09628 


. 1033 

.01201 

1.00262 

.1111 


.1033 

. 01402 

1.00337 

.1261 
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of all five fuel-filled spheres, the fuel was removed. Each sphere was then rinsed sev- 
eral times with demineralized water, vacuum dried, and refilled with the desired boric 
acid solution. An identical volume and reactivity worth determination was then made for 
each boric acid concentration as had been made for the fuel. 

The void reactivity worth correction for the difference in bubble size between fuel- 
and boric -acid-filled spheres was determined by use of the ZPR-I core center void worth 
data of appendix C. This correction, as a percent of the total boron -for -fuel substitution 
reactivity worth, amounted to less than 0. 5 percent. 

The following sequence of measurements is representative of the procedure used in 
determining the boron-for-fuel substitution reactivity worth, for a five-ball series, in a 
ZPR-I base core: 

(1) Fuel solution is pumped into the reactor tank until criticality is attained. Critical 
solution height and solution temperature are recorded. Solution critical height is cor- 
rected to 20° C by the temperature coefficient of appendix D. 

(2) A fuel-filled plastic ball is held at the core center by the holder of figure 1. Fuel 
is pumped into the reactor tank until the core is critical. Solution critical height and tem- 
perature are recorded. The critical height is corrected to 20° C. (This step is identical 
for each of the five balls. ) 

(3) Step 1 is repeated after the experiments indicated in step 2 are completed. The 
difference between the critical heights of steps 1 and 3 is due to evaporation of water 
vapor from the fuel solution during the time interval of the series of measurements. The 
critical height measurements of step 2 are corrected accordingly, if evaporation is as- 
sumed to be a linear function of time. The rate of change in solution critical height, due 
to water evaporation, ranged from about 0. 001 centimeter per hour for the 24. 892- 
centimeter base core to approximately 0. 005 centimeter per hour for the 67. 310- 
centimeter base core. Accordingly, the correction was small for these experiments. 

Critical heights of a critical core with the boric acid solution substituted for fuel so- 
lution were determined in an identical manner as outlined previously. The difference be- 
tween the critical height of a fuel-filled and a boric-acid-filled ball for each ball size and 
boric acid concentration is converted to reactivity worth by the method outlined in appen- 
dix B. 


Specific Reactivity and Base Core Flux Integral 

Since equation (12) is in a form suitable for evaluation, its components may now be 
considered. The reactivity worth of a boron-for-fuel substitution for each of the three 
boric acid concentrations in the 45. 720 -centimeter base core is plotted as a function of 
substitution volume in figure 2. For small volume insertions, reactivity is a linear func- 
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0 4 8 12 16 20 24 28 32 

Volume of boric acid solution, cm 5 


Figure 2. - Reactivity worth of unit spherical volume of boric acid solution of 45. 720- 
centimeter critical ZPR-I core at 20° C. Actual delayed neutron fraction for 
uranium 235, 0.0064. 

tion of substitution volume. The slope of the function, specific reactivity p_ (in units of 

o S 

(AK/K)/cm based on 0 = 0. 0064), is drawn for the three boric acid concentrations. 

The p thus determined is then plotted as a function of boric acid S A (0. 0253 eV) in fig- 
s a 

ure 3. The S^(0. 0253 eV) of the fuel solution of this base core is 0. 1191 per centimeter. 
A matching boric acid S A (0. 0253 eV) (see fig. 3) corresponds to p g of 0. 0001237 per 
cubic centimeter. Specific reactivities of boron -for -fuel substitutions performed in the 
other base cores, which were determined in the same manner as stated previously, are 
given in table V. 

In order to evaluate the flux integral of equation (12), the flux <h(r) and adjoint flux 
<t>*(r) distribution functions must be known. This implies that the extrapolation distance 
6 to zero neutron flux beyond the core boundaries is also known. Activation of thin 

12 


j 




l/l 

Q_ 

fc 

OJ 


o 



138x10 

134 


-.130, 


-. 126 


122 


-.118 


-.114 


-. 110 1 

.1020 



Figure 3. - Specific reactivity of boric acid solutions in center of 45.720-centimeter critical core, as function of 
boric acid solution macroscopic neutron absorption cross section. 


TABLE V. - EXPERIMENTAL RESULTS 21 


Base 

Uranyl 

Specific reactivity 

Absorption 

Core flux 

Thermal 

Core flux 

Effective 

Effective 

core 

fluoride - 

measurements 

cross 

integral 

power 

integral 

delayed 

to actual 

critical 

water 



section 

(dysprosium), 

fraction, 

(total), 

neutron 

delayed 

height 

solution 

p S ' 

Q 

P s /6, 
. «2 

mismatch 

f 3>(r)<$*(r)dr, 

f 

f $(r)<$*(r)dr, 

fraction, 

neutron 

at 20° C, 

hydrogen to 

(AK/K)/cm J 

$/cirr 

correction, 

J v - — 


y V “ - “ 

^eff 

ratio, 

cm 

uranium 235 



6A 

cm 3 


cm 3 


y 


atom ratio. 










R 


(b) 





(b) 

(b) 

24. 892 

190. 1*0. 6 

-0. 0002009 

-0.03139*0.00012 

0.00139 

3221±229 

0.910 

3540±251 

0. 0090±0. 0006 

1.40*0. 10 

45.720 

473. 1*3. 6 

-.0001237 

-. 01933*0. 00008 

-.00057 

5802±342 

.960 

6043±357 

. 0086*0.0005 

1.34*0.08 

67.310 

565. 8*4. 8 

-.00008540 

-. 01334*0. 00005 

-.00121 

8839±468 

.965 

9159±485 

.0082*0.0004 

1. 28*0.07 


a All error limits are 95 percent confidence limits. 
b 0 = 0. 0064 = 1 $. 


uranium 235 detectors within a base core provided the correct <&(r) distribution function 
and core boundary extrapolation distance. The uranium detectors used were 0. 635- 
centimeter-diameter and 0. 0279-centimeter-thick uranium -aluminum alloy foils, which 
contained 15. 5 percent uranium by weight. The uranium is 93. 3 percent enriched in the 
uranium 235 isotope. Only one core traverse with uranium 235 detectors was necessary 
since axial and radial mapping in all of the base core sizes had previously been done with 
dysprosium detector foils. It will be shown that the flux integral obtained from the dys- 
prosium foil activity distribution can be corrected to the flux integral obtained from the 
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Relative neutron flux distribution, <£(z) 


1.2 


1.0 


■ 4 


.2 



I I I I I I I 
Alloy foils 

0.635-cm-diam uranium 235-aluminum 
0. 318-cm-diam dyprosium-aluminum 

Relative neutron flux distribution, 
<£(z) 

cos(0. 27819 z H Determined by least- 
cos{0. 28964 z) J squares fit of data — 

I I 

Extrapolation^ 

distance, 

6 , 

cm 

1. 87±0. 19- 
1. 31±0.19 


-4 -2 0 2 4 6 8 

Axial distance from core center, z, cm 
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Figure 4. - Determination of extrapolation distance of neutron flux distribution along vertical axis of ZPR-I core 
using uranium 235 and dysprosium foil activations. Base core height, 24.892 centimeters. 



Figure 5. - Extrapolation distance beyond axial and radial core boundaries as function of ZPR-I critical core 
solution height. Error limits, ±1 standard deviation. 
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uranium 235 distribution. The dysprosium detectors used were 0. 3 175 -centimeter- 
diameter and 0. 0127-centimeter-thick dysprosium -aluminum alloy foils. The alloy was 
5 percent dysprosium by weight. 

Foil activation data were fitted in the axial and radial directions by cos(ttz/Z) and 
J Q [2. 405(r/Y)] functions, respectively. Both uranium 235 and dysprosium axial distribu- 
tions for the 24. 89 2 -centimeter core are plotted in figure 4. A least-squares fit of the 
two sets of data by a cos(irz/Z) function resulted in a value of 6 of 1. 31 centimeters by 
the dysprosium foil activations and 1. 87 centimeters by the uranium 235 foil activations. 
Similar neutron flux distributions for the other base core sizes were obtained. A plot of 
the resulting 6 as a function of critical core height is shown in figure 5. 

K it is assumed that 6 is energy independent and $(r) and <f>*(r) of equation (12) 
have the same spatial distribution, the product <&(r)$*(r) would then be equal to <f>^, 

which is the square of the experimentally measured neutron flux distribution. To verify 
2 

a dependence of reactivity worth variation over the core volume, a boric -acid -filled 
sphere, having the same £^(0.0253 eV) as the core fuel solution, was traversed both ax- 
ially and radially in the 24. 892-centimeter core. The largest sphere (28. 30 cm^ volume) 



0 2 4 6 8 10 12 14 

Axial distance from core center z cm 


Figure 6. - Relative neutron flux distribution as function of distance along 
vertical axis from center of a critical cylindrical core. Base core height, 
24.892 centimeters. 


was used in order to provide mea- 
sureable reactivities of sufficient 
magnitude in core regions of least 
importance. However, the combi- 
nation of large sphere size and 
small core volume enhanced the 
positional errors. The resulting 
relative reactivity data points, with 
their estimated positional errors, 

are plotted in figures 6 and 7. The 

2 

data appear to follow the $ dis- 
tribution reasonably well. The 
flux integral of equation (12) was 
evaluated for the base core sizes, 
utilizing the measured 6 and flux 
distributions discussed previously. 

The flux distributions obtained 
from dysprosium foil activations 
were essentially thermal since the 
cadmium ratios (ref. 17) of the ac- 
tivated dysprosium foils in the base 
cores were greater than 100. 
Therefore, the uranium 235 foils 
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Figure 7. - Relative neutron flux distribution as function of distance along radius 
from center of critical cylindrical core. Base core height, 24. 892 centimeters. 


indicated the total fission rate while dysprosium foils indicated only the thermal fission 
rate. The ratio of the flux integral using the dysprosium extrapolation distance to the 
flux integral using the uranium 235 extrapolation distance should then be representative 
of the ratio of the thermal fission rate to the total fission rate f . This was verified to 
within 3 percent with a diffusion theory calculation assuming a dysprosium "cutoff" of 
about 3 electron volts. The uranium 235 traverse was performed only for the 24. 892- 
centimeter core and calculated values of f were then used for the other two cores. The 
values of f and the total flux integral are given in table V (p. 13). These flux integrals 
have been normalized according to equation (11). 


Experimental Determination of Effective Delayed Neutron Fraction 

The results presented in table V include the measured specific reactivities of the 
boron-for-fuel substitutions performed in the three base cores and the corresponding flux 
integrals. The flux integrals that utilized dysprosium foil activations for relative neutron 
flux distributions represent the contribution of thermal neutrons only. Division of this 
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Figure 8. - Effective to actual delayed neutron ratio as function of hydrogen to uranium 235 atom ratio 
(95 percent confidence limits). Actual delayed neutron fraction for uranium 235, 0.0064. 


integral by the thermal power fraction f results in the total core flux integral. This 
total flux integral weights effectively the specific reactivities by a power distribution that 
includes neutrons of all energies. The effective delayed neutron fraction is calcu- 

lated from equation (12) using the specific reactivities and total flux integrals for the base 
cores of table V. Shown in table V is the experimentally determined ratio y = 'W' 3 £or 
the three base cores. Values of y, as a function of R, are also plotted in figure 8. 

The accuracy of the determination of depends on the precision of the specific 

reactivity p_ and the base core flux integral measurements. The standard deviation s 
of p_ is less than ±0. 2 percent for all boron-for-fuel substitutions. The standard devia- 
tion of the measurement of critical core extrapolation distance 6 is ±0. 19 centimeter. 

A 95-percent confidence limit of ±0. 38 centimeter (2s) applied to the extrapolation dis- 
tances results in error limits of ±29, ±22, and ±12 percent for the ZPR-I reactors having 
24.892-, 45.720-, and 67. 310-centimeter heights, respectively. The corresponding er- 
ror limits in the base core flux integrals are ±7. 1, ±5. 9, and ±5. 3 percent. The error 
limits of |3 e jj are, therefore, essentially the same as the error limits of the base core 
flux integrals. The standard deviations in the determination of the hydrogen to uranium 
235 atom ratio R for the three base cores 190. 1, 473. 1, and 565. 8 are ±0. 16, ±0. 38, 
and ±0. 42 percent, respectively. 


THEORETICAL CALCULATIONS FOR EFFECTIVE DELAYED NEUTRON FRACTION 

Theoretical calculations for are made for each reactor using the ’’parallel 

group” technique described in reference 7. This technique is particularly useful when 
multigroup spatial computer programs are available. Usually the number of neutron en- 
ergy groups required is almost double the number normally required for a calculation of 
the effective multiplication factor. For bare homogeneous reactors some simplification 
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table vi. - energy group structure is possible with little loss in accuracy. 

The criticality calculations are performed 
using five neutron energy groups. The linear 
nature of the equations allows the contributions 
from the delayed neutrons to be written as sepa- 
rate equations and solved simultaneously, in 
parallel, with the prompt equations. In this par- 
ticular group structure, given in table VI, no 
delayed neutrons appear in group 1; hence, nine 
equations result. Reference 7 reports the use 
of four prompt groups and four delayed groups 
with only one energy group transfer. The extension to more groups with a complete down- 
scattering transfer matrix is straightforward. These equations, which are written in dif- 
fusion theory form, are the following: 


Group 

Energy 

Uranium 235 

Prompt 

fission 

fraction 

Delayed 

fission 

fraction 

1 

14.92 to 0.498 MeV 

0.868659 

0 

2 

0.498 to 0.183 MeV 

.098832 

1.0 

3 

183 to 11.7 keV 

.032508 

0 

4 

11 700 to 0.414 eV 

0 

0 

5 

0.414 to 0 eV 

0 

0 


-dPv 2 *p + ( S P + sP _ 2 + sP _ 3 4 sP _ 4 4 sP. 5 )*P = -A s 


x p 

_3 

K 


- D Pv 2 *P 4 (sP 4 sP_ 5 )*P = S P_ 4 ®P 4 sP_ 4 *P 4 S P_ 4 *P 
-DjV 2 <H> 4 2 P 4 P = sP_ 5 *P 4 S P_ 5 *P 4 sP_ 5 *P 4 sP_ 5 *P 


.d„2 (T> d 

) 7 v <6, 


(4 


,d 

'7—8 




■? 


■pvdyy 2 id / •yvd \ id _ d id yd 1 1 

" D 8 V ^8 + l s 8 + ^8—9/8 " S 6-8 $ 6 + S 7-8^' 



(13a) 

fs4 S P_ 2 *P 

(13b) 


(13c) 

4*2 + S 3-4*3 

(13d) 

5*S - 

(13e) 

= — S 
6 K 

(13f) 

r>d id 

S 6-7^6 

(13g) 

yd id 

' S 7-8^7 

(13h) 
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2 t d ydaid yd yd t d yd id 

-°9 V *9 + S 9*9 = S 6-9*6 + S 7-9*7 + Z 8-9*8 


(13i) 


where 

5 9 

S = Y vZf ®P + V vzf 
jL/ f, g g Zj f, g g 
g=l g=6 

and 


These nine equations are solved using a one -dimensional diffusion theory program 
for the IBM 7094 computer system. The macroscopic cross sections for the first four 
groups are averages over a spectrum resulting from the prompt spectrum of uranium 235 
fission neutrons as obtained from a proprietary program GAM-II (ref. 11). The 
TEMPEST program (ref. 12) is used to obtain the macroscopic cross sections for the 
thermal groups 5 and 9. A Wigner -Wilkins spectrum is assumed. The cross sections 
for the delayed groups 6, 7, and 8 are obtained from GAM-II using a fission spectrum of 
delayed neutrons with a mean energy of 450 keV (refs. 13 and 14). 

The first solution of the nine equations is obtained using the following fission fraction 
values: 


x i = x i 
x P = X 2 -/3 

X 3 = x 3 

x 6 = ^ 

where y = 1.0 and the are given in table VI. The resulting eigenvalue is K^. 

Next, the same set of equations is solved with a perturbation to Xgj that is, Xg = 2/3 and 
y Xj = 1 + /3. The eigenvalue is Kg, and /3 pf ^ //3 is given by 


fleet _ K 2 - K 1 

0 KgKj/3 
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For bare homogeneous reactors, it is not necessary to use the preceding method. 
Since the delayed neutrons have the same spatial distribution as the prompt neutrons and 
the spatial solutions are cosine and Bessel functions, the ratio /3 eff //3 can be obtained by 
solving only the last four equations with Xg * 1. The eigenvalue maybe interpreted as 
j3 pff //3. The critical geometry is assumed such that a solution of the first five equations 
gives K = 1. To be consistent with the previous definition of 0 e ff/0, the following rela- 
tions are used: 


^eff _ K 2 ~ K 1 

0 KgK^ 


for nine group cases, and 


Paff _ <*¥ + K l> ~ K 1 

0 (K 2 0 + K^Kj/3 


(K 2 0 + KjJKj 


for four group cases. 

A summary of the calculations is given in table VII. Case 1 is based on the nine 

TABLE VH. - CALCULATED RESULTS FOR EFFECTIVE TO 
ACTUAL DELAYED NEUTRON RATIO 


Case 

Ratio of 
hydrogen and 
uranium 235 
atom density, 
R 

Effective multiplication 
factor of a reactor 

Effective to actual delayed 
neutron ratio, 

Peff/ p 

First, 

K i 

Second, 

K 2 

Theoretical 

Experimental 

1 

190 

1 . 000401 

1.008999 

1.3268 

1.40 


473 

1.001275 

1.009484 

1.2650 

1.34 


562 

1.000821 

1.008900 

1.2463 

1.28 

2a 

190 

1.000169 

1.3413 

1.3294 

1.40 


473 

1.001085 

1.2809 

1.2677 

1.34 


562 

1.000647 

1.2606 

1.2489 

1.28 

2b 

190 

1.000679 

1.3946 

1.3804 

1.40 


473 

1.000473 

1.3200 

1.3077 

1.34 


562 

1.000464 

1.2976 

1.2857 

1.28 
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r 


group solution for each of the three reactors. Case 2 is based on the four delayed groups 
only. Case 2a, as well as case 1, uses energy-dependent extrapolation distances 
6(E) = 0. 71 Xj r , while case 2b uses an equivalent energy independent extrapolation dis- 
tance 6 which is equivalent in the sense that 6 is chosen so that criticality is computed 
to approximately the same accuracy as when using 6(E). 

Cases 1 and 2a are in agreement and both are calculated with an energy-dependent 
extrapolation distance. Case 2b, if energy -independent extrapolation distances are used, 
gives significantly higher values of ^ e ff//3. It is this case that is most directly compara- 
ble with the experiment, since energy independent extrapolation distances are assumed in 
the experiment also. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 7, 1966, 

122-28-03-02-22. 
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APPENDIX A 


c 

D 

E,E' 

f 

m,\ 

f(E’)J 

H 

AH 

h 

i 

J 0 

j 

K 

K V K 2 

AK/K 

i 

N 

P(r) 

P*(r) 


R 

r 

r 

S 


SYMBOLS 


transverse buckling, cm 
constant 

diffusion coefficient, cm 
energy, eV 

thermal power fraction ratio of thermal fission rate to total fission rate 
fraction of fission neutrons per unit energy 


modified critical height of core, h plus 9. 8 cm, cm 
incremental change in critical core height, cm 
actual critical height of a ZPR-I core, cm 
i in substitution 

Bessel function of first kind of zero order 
j group 

effective multiplication factor of a reactor 

eigenvalues 

(K - 1)/K 

mean prompt neutron lifetime, sec 

0 

atom density, atom/cm 


neutron production rate at position r per cm"*, (r)<f>(r) for g energy 

groups of neutrons g & 

fission spectrum weighted adjoint at position r per cm'*, ^ x g $|(£) for g en " 
ergy groups of neutrons g 

ratio of hydrogen to uranium 235 atom densities 

distance along a radius perpendicular to axial core centerline, cm 

position vector 

fission source including both prompt and delayed neutrons 
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s standard deviation 

T stable reactor period, sec 

it. 

Tj stable reactor period corresponding to i substitution, sec 
At incremental change in fuel solution temperature 

V total core volume 

X neutron fission fraction for given energy group 

Y actual critical core radius plus 6 

Z actual critical height plus 26 

z distance from core center to point on cylinder axis, cm 

6 

/3 actual delayed neutron fraction for uranium 235, ^ /3^ 

j=l 

/3 g ff effective delayed neutron fraction (see eq. (12)) 

iy. 

/3. fraction of total fission neutrons in j delayed group 
y effective to actual delayed neutron ratio, /3 e ff/ 0 

6 energy independent extrapolation distance from reactor core boundary to point 
at which flux is zero 

6A absorption cross section reactivity mismatch correction (see eq. (7)) 

6(E) energy dependent extrapolation distance, cm 

6 p T incremental reactivity worth of substitution of nonfissionable material for 

1 fissionable material in core region of critical reactor 

62 incremental perturbation in 2 

a a 

Xj decay constant for j th group of delayed neutrons, sec -1 
X tr transport mean free path, cm 

i>2j number of neutrons per fission times fission cross section, neutrons/cm 
p reactivity, AK/K based on /3 = 0. 0064 

A p incremental change in reactivity associated with incremental change in critical 

core height, AH 

3 

p_ specific reactivity, reactivity per cm of boron-for-fuel substitute 
s 
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S A (0. 0253 eV) macroscopic neutron absorption cross section at 0. 0253 eV, cm 


-1 


S g-~g+n 

a A (0. 0253 eV) 

T i 

*2 

$(r) 

**(r) 

<&(r, E) 

<£>*(r, E) andl 
<&% E’) J 

Subscripts: 

calc 

exp 

S 

Superscripts: 

d 

P 


macroscopic absorption cross section for a given energy group, cm 

transfer cross section from energy group g to group g + n 

(n = 1, 2, . . . and g + n =£ total number of energy groups), cm' 

microscopic neutron absorption cross section at 0. 0253 eV, b 

volume region of i in substitution 

square of experimentally measured neutron flux distribution 
neutron flux at position r 
adjoint neutron flux at position r 
neutron flux per unit energy at r and E 

adjoint per unit energy at r and E or E’ 

calculated values 
experimentally determined values 
energy group numbers, 1, 2, . . . , 9 

delayed neutrons 
prompt neutrons 
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APPENDIX B 


CALIBRATION OF AN NASA ZPR-I CRITICAL CORE 
Introduction 

The ZPR-I calibration curve to be derived for an unperturbed core gives directly the 
reactivity worth of an increment of fuel solution height due to a uniform change of the 
uranyl fluoride - water fuel solution 2^(0. 0253 eV). This same calibration curve can be 
used to determine the reactivity worth of an incremental change in fuel height due to a 
small core-centered perturbation, up to a certain sample size and reactivity worth. Sam- 
ple sizes up to nearly 4 centimeters in diameter that have a reactivity worth of nearly a 
dollar can still be related to reactivity by the unperturbed core calibration curve. The 
usefulness of these findings are twofold. First, absolute reactivity worths of moderate 
core centered perturbations (reactivity worth less than 1 $) can be found easily in a 
ZPR-I core. Second, the reactivity worth of a perturbation in a certain base core con- 
figuration can be found although the measurement may have been performed in a critical 
core of a slightly different height. 


Slope of Unreflected Unperturbed Critical Uranyl Fluoride - Water Solution Core 

Reactivity p can be calculated from a measured stable reactor period T by the in- 
hour equation (ref. 15) 


6 

P = - + rY — ^ — (Bl) 

T L~i 1 + X.T 
j=l 3 

The stable reactor period associated with an incremental increase of fuel solution 
above a critical core height has been measured for several fuel concentrations in ZPR-I. 
Small increments of fuel were added to a critical reactor to obtain stable reactor periods 
that were always greater than 200 seconds. For the ZPR-I reactors considered, H 
varied from about 25 to 60 microseconds. As a result, the reactivities determined can 
be considered independent of mean prompt neutron lifetime £. Moore’s reactivity table 
(ref. 16), derived from the inhour equation with /3 = 0. 0064 and y= 1.0, was used for 
convenience in determining reactivities from the measured stable periods. An example 
of such a determination is shown in figure 9. The slope of this curve is Ap/AH. 
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25.060 25.070 25.080 25.090 25.100 25.110 
Solution height, cm 


Figure 9. - Reactivity as function of uranyl 
floride - water fuel solution height above 
unperturbed critical core. Actual delayed 
neutron fraction for uranium 235, 0.0064. 


When the uranyl fluoride - water solution critical heights are plotted as a function of 
(Ap/AH) on linear graph paper, the data can be fitted by a straight line that inter- 
cepts the ordinate a few centimeters below zero critical height. Gwin (ref. 5) presents 
Ap/AH data in this manner for several different configurations of the uranyl fluoride - 
water critical solution geometries. His data extrapolated to an ordinate intercept of 
-7. 5 centimeters. A fit of ZPR-I Ap/AH data, including Gwin’s data from reference 5 
for a 30. 5-centimeter-diameter critical uranyl fluoride - water core, resulted in an ordi- 
nate intercept at -9. 8 centimeters as shown in figure 10. If h equals the actual ZPR-I 
core critical height, then let 


H = h + 9. 8 cm 


(B2) 


so that 
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Figure 10. - Critical height as function of (reactivity per unit change in core 

height)' 1 ^, (Ap/AH)" 1 ^. Actual delayed neutron fraction for uranium 235, 
0. 0064. 


H = 



(B3) 


The constant C in equation (B3) is evaluated from figure 10 and equals 7.4303. 
Equation (B3) can be written as 


Aa. = 410. 22 H" 3 (B4) 

AH 


In the limit as A p and AH approach zero, equation (B4) can be written as the dif- 
ferential dp, where 
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(B5) 


dp = 410.22 H" 3 dH 


which can be integrated to obtain reactivity. 


Slope of Unreflected Perturbed Critical ZPR-I Cores 

A ZPR-I reactor of one critical height can easily be changed to a core of a different 
critical height by simply adding water to or evaporating water from the uranyl fluoride - 
water solution. This action changes the hydrogen to uranium 235 atom ratio R and 
S A (0. 0253 eV) of the fuel solution. The Ap/AH slope of any particular core height can 
be found from figure 10. If the fuel solution were altered slightly by the uniform addition 



(Ap/AH )' 1 ' 3 

Figure 11. - Critical height as function of (reactivity per unit change in core 
height)” 1 ^ 3 , (Ap/AH)“ 1/3 , for unreflected, perturbed, cylindrical reactors of 
uranyl fluoride-water solution in aluminum containers. Actual delayed 
neutron fraction for uranium 235, 0.0064. 
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of a neutron absorber, for example, boric acid to the fuel solution, and the Ap/AH slope 
of this perturbed reactor measured, then this slope is expected to be identical to that of 
unperturbed core of the same critical height but of a different R ratio. 

A comparison of a localized absorber with a uniform dispersion of absorber is de- 
sired. A Ap/AH slope was measured for several ZPR-I critical cores with core- 
centered perturbations. The results are shown in figure 11, which is basically a repro- 
duction of figure 10 but with the data points of core-centered perturbations superimposed 
on the unperturbed critical core curve. With the exception of the large cadmium box ab- 
sorber, the core-centered absorbers shown have Ap/AH slopes nearly identical to un- 
perturbed critical cores of the same height but a different uranyl fluoride - water solution 
concentration. The two 28-cubic-centimeter spheres of boric acid solution, of the same 
S A (0. 0253 eV) as the uranyl fluoride - water solutions they were inserted into, represent 
perturbations at the core center with reactivity worths as high as 0. 006 AK/K. The seven 
2. 54-centimeter-diameter foils of boron 10 (S^(0. 0253 eV) = 1. 38) provided a rather 
severe core-centered thermal neutron flux perturbation with a reactivity worth of approxi- 
mately 0. 0038 AK/K, yet the measured Ap/AH slope was identical to that of an unper- 
turbed core of the same height. The size of the cadmium box (5- by 3. 8- by 0. 64-cm) 
and reactivity worth (about 0. 0089 AK/K) apparently distorted the thermal neutron flux 
too severely to be classified as a simple core perturbation. 

The specific reactivity p_ of the boron -for -fuel core-centered substitution per- 
formed for this report depended strongly on the reactivity worth measurement of the 

Q 

smallest ball used (approximately 4. 6 cm ) for determination of unit volume reactivity 
worth. The measured reactivity worth of this boric-acid-filled ball was much less than 
0. 001 AK/K for all ZPR-I base core sizes. Based on the data presented in figure 11, it 
can be seen that use of the H against (Ap/AH) functional relation to obtain reactivi- 
ties is valid for perturbations with reactivities up to approximately one dollar. 

Integral Calibration Curve Relating Incremental Change in 
Core Height to Reactivity 

The foregoing discussion and data of figure 11 have shown that the reactivity associ- 
ated with an incremental change in ZPR-I critical core height AH caused by a uniform 
change in the uranyl fluoride - water solution concentration, is equivalent to a AH 
caused by a moderate core centered perturbation. Therefore, a series of core perturba- 
tions made in cores that differ slightly in critical height from day to day, because of the 
continuous evaporation of water from the fuel solution, can be normalized to a particular 
base core critical height. An integral calibration curve to normalize the unperturbed 
critical core heights, and to determine the reactivity worth of a core perturbation, can 
be obtained by integration of equation (B5): 
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(B6) 


p = /dp = 410. 22 /h - 3 dH + C 

As an example, consider the 45.720-centimeter-base critical core. Integration of 
equation (B6) gives 


p = -(205. 11)H" 2 + C (B7) 

The constant C is evaluated by setting p = 0 at H = 45. 720 + 9. 8 = 55. 520 centimeters 
so that 


C = - 205 ' 11 = 0. 066541 (B8) 

(55.520) 2 

and the equation for p as a function of H is 

p = -(205. 11)H“ 2 + 0. 066541 (B9) 

As an example of the use of equation (B9) for the 45. 720-centimeter-base core 
(H = 55. 520 cm), one reactivity determination of a boron -for -fuel substitution will be 
given. The unperturbed core critical height at 20° C is 45. 8384 centimeters 
(H = 55. 6384 cm). The reactivity worth of the increment of fuel solution above the base 
core height is 0. 000282 AK/K. Next, the critical height at 20° C with a core centered 
4. 5296-cubic-centimeter ball of boric acid solution (2^(0. 0253 eV) = 0. 1229 cm - *) is 
46. 3098 centimeters (H = 56. 1098 cm). The reactivity worth of fuel solution and core 
perturbation is equal to 0. 001392 AK/K. Thus, the net reactivity worth of the substitu- 
tion equals 0. 001392 - 0. 000282 = 0. 001110. 

Reactivities associated with core centered perturbations in the 24. 892- and 67. 310- 
centimeter base critical cores were found in the same manner as for the 45. 720- 
centimeter core. The equations for these two cores are 

p = -(205. 11)H~ 2 + 0. 170424 (B10) 

for the 24. 892-centimeter core and 

p = -(205. 11)H" 2 + 0.034496 (Bll) 

for the 67. 310-centimeter core. 
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APPENDIX C 


REACTIVITY WORTH OF ZPR-I CORE CENTERED VOID 

As described in the EXPERIMENTAL PROCEDURE section, a small air bubble re- 
mained in each of the boric-acid- and fuel-filled plastic spheres of table n (p. 8). There- 
fore, it was necessary to measure the reactivity worth of a core -centered void volume in 
order to obtain a correction for this difference in bubble size. This measurement was 
accomplished by determining the core-centered reactivity worth of each of the empty 
plastic spheres. An example of one base core measurement of this void worth is shown 
in figure 12 for the 45. 720 -centimeter core. The reactivity per unit volume for small 

_C 

volumes is a linear function of volume with a slope of 5. 21X10 . Figure 13 shows the 

central core void volume worth of all three base cores. 



Figure 12. - Worth of void in center of 45.720-centimeter (18.000 in. ) critical ZPR-I core 
at 20° C. Actual delayed neutron fraction for uranium 235. 0.0064. 


31 



26 30 34 38 42 46 50 54 58 62 66 

Critical height, cm 

Figure 13. - Core-centered unit volume void worth as function of ZPR-I critical core height at 20° C. 
Actual delayed neutron fraction for uranium 235, 0.0064. 



APPENDIX D 


TEMPERATURE COEFFICIENT OF REACTIVITY FOR AN NASA ZPR-I REACTOR 

The ZPR-I is operated at very low power and therefore undergoes no temperature 
change as a result of heat generation within the core. The variation of fuel temperature 
is caused by the change in the reactor room air temperature. The reactor room is pro- 
vided with a heating system and an exhaust fan for ventilation but not with any air condi- 
tioning equipment. Day to day changes in room air temperature are usually slight. 

Longer term seasonal temperature variation causes gradual room air temperature change 
from approximately 20° to 25° C. 

For a reactor temperature coefficient determination, the critical height of the ZPR-I 
fuel solution at normal room temperature is first measured. The reactor room air is 
then heated until approximately a 2° to 5° C increase in fuel solution temperature occurs. 
The new critical height of the warmer fuel solution is then measured. The ratio of the 
incremental difference between these two critical heights AH to the number of degrees 
change in solution temperature At that had occurred is the temperature coefficient. 

This coefficient has been found to be a constant over the temperature range involved 
within the accuracy of the measurement. 

The temperature coefficient AH/ At was measured for three base cores having 
critical heights of 24.892, 45.720, and 67.310 centimeters. These temperature coeffi- 
cients as functions of ZPR-I critical fuel solution height are plotted in figure 14. Tem- 
perature coefficients in reactivity per degree centigrade Ap/At for the same base cores 
were obtained by converting the AH of figure 14 to A p from figure 10 (p. 27). The re- 
sulting reactivity temperature coefficient as a function of critical height is plotted in fig- 
ure 15. 
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Temperature 

coefficient, Ap/At Temperature coefficient, AH/ At, cm/°C 




24 26 30 34 38 42 46 50 54 58 62 66 70 

Critical height, cm 

Figure 15. - Temperature coefficient obtained from figures 10 and 14. Actual delayed neutron fraction for 
uranium 235, 0.0064; temperature range, 20° to 26° C. 
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